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Dark matter that was once in thermal equilibrium with the Standard Model is generally prohibited 
from obtaining all of its mass from the electroweak or QCD phase transitions. This implies a new 
scale of physics and mediator particles needed to facilitate dark matter annihilations. In this work, 
we consider scenarios where thermal dark matter annihilates via scalar mediators that are colored 
and/or electrically charged. We show how partial wave unitarity places upper bounds on the masses 
and couplings on both the dark matter and mediators. To do this, we employ effective field theories 
with dark matter as well as three flavors of sleptons or squarks with minimum flavor violation. For 
Dirac (Majorana) dark matter that annihilates via mediators charged as left-handed sleptons, we find 
an upper bound around 45 TeV (7 TeV) for both the mediator and dark matter masses, respectively. 
These bounds vary as the square root of the number of colors times the number of flavors involved. 
Therefore the bounds diminish by root two for right handed selectrons. The bounds increase by root 
three and root six for right and left handed squarks, respectively. Finally, because of the interest in 
natural models, we also focus on an effective held theory with only stops. We find an upper bound 
around 32 TeV (5 TeV) for both the Dirac (Majorana) dark matter and stop masses. In comparison 
to traditional naturalness arguments, the stop bound gives a firmer, alternative expectation on when 
new physics will appear. Similar to naturalness, all of the bounds quoted above are valid outside of 
a defined hne-tuned regions where the dark matter can co-annihilate. The bounds in this region of 
parameter space can exceed the well-known bounds from Griest and Kamionkowski [I]. We briefly 
describe the impact on planned and existing direct detection experiments and colliders. 


I. INTRODUCTION 

Understanding the nature of dark matter is one of the most pressing, unresolved problems in particle physics. Dark 
matter is needed to understand structure formation, the observed galactic rotation curves m and the acoustic peaks 
in the cosmic microwave background [5]. Moreover, the dark matter relic abundance is measured to be [5] 

h 2 Q c = 0.1199 ±0.0027. (1) 

A compelling argument for the origin of this abundance is to assume dark matter was once in thermal contact with 
the baryon-photon plasma during the early universe. Since all known forms of matter in the universe were once 
in thermal equilibrium, this type of dark matter is theoretically persuasive. In this scenario, the measured relic 
abundance is controlled by dark matter annihilations into Standard Model (SM) particles. Because of constraints 
from the observed large scale structure in the universe, dark matter must be stable and non-relativistic (cold) when 
departing thermal equilibrium J3]. 

The Standard Model (SM) alone cannot account for the missing matter in the universe [6]. Current experi¬ 
mental constraints, however, provide some guidance on the structure of the underlying theory. For example, the 
lack of large missing energy signatures at the Large Hadron Collider (LHC) [7HT7] and other colliders p~8l [25] 
suggests that the dark matter is either heavy or has very small couplings with the SM so that it is not produced in 
high-energy collisions. Additionally, direct detection experiments [261428] . updated precision electroweak constraints, 
and precision Z-pole experiments [291I3T] all severely constrain the direct coupling of dark matter to the SM Higgs 
and/or Z bosons. These constraints all imply that dark matter cannot obtain all of its mass from the SM Higgs 
alone [30] . Mediator-facilitated interactions help to evade current experimental constraints by partially decoupling 
the dark matter from the SM. Should these scenarios be realized in nature, the discovery of the mediator particles 
would be an important step in understanding the nature of dark matter. It is therefore crucial to place bounds 
on the masses and couplings of these mediators. The most popular ways for dark matter to annihilate via a 
mediator particle are through scalars that are colored or charged, the Higgs boson [32] . or via a new neutral gauge 
boson. Some of us explored the perturbative unitarity constraints on the Higgs [331435] and the gauge portals [36] . 
In this work we focus on constraining scenarios where fermionic dark matter annihilates via charged or colored scalars. 

In supersymmetric (SUSY) models, the charged and colored scalars are sleptons and squarks, respectively. 
We employ an effective field theory including these particles as a benchmark in our analysis. In addition, because of 


2 


the intense focus on natural SUSY models, we particularly consider an effective field theory with dark matter and the 
supersymmetric partner of the top quark (stop). The stop is the key for achieving natural, supersymmetric models of 
new physics. Note naturalness [37] has been the dominant paradigm that has shaped models of new physics for the 
last 35 years. It is minimally satisfied with a relatively light stop [381 which cancels the quadratic divergences to the 
SM Higgs mass from the top quark. Moreover if the neutralino dark matter is only one component of the measured 
relic abundance (e.g., axionic plus neutralino dark matter), then the constraint on the stop mass improves. Thus, 
our bounds may be useful in shaping viable R-parity conserving supersymmetric models. 

We apply unitarity constraints in a manner reminiscent of Griest and Kamiokowski pQ. However, there are 
important differences: here we focus on perturbative unitarity constraints which determine, in particular, when the 
dark matter couplings become strong. Importantly, weakly interacting, massive dark matter and perturbativity have 
always had an intertwined conceptual association. In models that violate perturbative unitarity, the elementary dark 
matter candidate is efficiently forming bound states as well as annihilating as the temperature decreases towards 
the thermal decoupling temperature. The annihilation diagrams for the elementary dark matter candidate can be 
easily altered (and often times dressed with lighter states) to form diagrams in which the bound state decays. The 
bound state will therefore have a lifetime much shorter than the age of the universe. Thus, in models that violate 
perturbative unitarity, the elementary dark matter candidate is not an asymptotic state and the bound state that it 
forms is not a viable dark matter candidate. Note, it is well known that viable dark matter candidates exist that are 
the result of strongly coupled or confining hidden sectors. However, in these models the dark matter annihilation 
processes are still perturbative [39] • We show our perturbative unitarity constraints are improved in comparison to 
the updated Griest and Kamiokowski bounds m- Of central importance is the fact that our methodology places 
constraints on any particle associated with the dark matter annihilation. Our bounds on the masses and couplings 
of the charged and colored scalars are novel. 

Our basic perturbative unitarity arguments are straightforward. Dark matter annihilating via the electrically 
charged and/or colored scalars does so via t-channel diagrams. As one increases the mass of the mediator, the 
annihilation cross section is diminished unless the Yukawa coupling (between the fermionic dark matter, SM fermions 
and mediators) is raised. A sufficiently large annihilation cross section is needed to satisfy the relic abundance 
constraint. This process also forces the dark matter mass to be raised as well. Eventually, the mediator mass is 
sufficiently large that the couplings violate perturbative unitarity therefore setting an upper bound on both the dark 
matter and mediator masses. The parameter space for this process is additionally constrained by requiring the dark 
matter relic abundance to be less than or equal to equation ([l]). Conversely, small mediator masses are constrained 
by collider searches and other measurements. 

In the next section, we introduce simple models onto which we place our unitary and relic abundance bounds. 
Section III. introduces our unitarity constraints. Section IV. computes the dark matter relic abundance and direct 
detection cross sections. Section V. describes existing experimental constraints from the LHC. We implement all 
of our constraints and perform a parameter scan in Section VI. The basic unitarity bounds appear here. Future 
experimental signatures are in Section VII. Conclusions and Appendices follow. 


II. REPRESENTATIVE SFERMION MODELS 

Throughout we consider effective theories where the only relevant degrees of freedom are a SM singlet fermionic 
dark matter particle (denoted by x and potentially having either a Dirac or Majorana mass), a scalar mediator, and 
the SM particles. For this class of models, we consider scalar mediators that facilitates the following tree-level DM 
annihilation processes, 


X + X-> f + f, (2) 

through a Yukawa-type coupling between the SM fermion /, the DM particle %, and the scalar mediator. As a result, 
/ must have the same SM quantum numbers as the scalar mediator. Because we consider all possible SM final states, 
we consider scenarios with all possible SM quantum numbers for the scalar mediators. In order to stabilize the DM, 
both it and the scalar mediator must be charged under the same dark matter stabilization symmetry. Because our 
scalar mediators are similar to the sfermions in SUSY models, we will therefore employ the SUSY nomenclature to 
describe them. However, the key difference with respect to SUSY theories is that we allow the DM-sfermion-fermion 
coupling to take an arbitrary value up to the perturbative limit. 
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We consider an effective scenario in which the DM-sfermion-fermion coupling is the only non-gauge interac¬ 
tion for the sfermion m- The sfermions will either be purely left-handed, /l, or right-handed states /r. Calling 
the DM-sfermion-fermion coupling A^, and considering a single set of sfermion quantum numbers, the interaction 
Lagrangian is 


£ — ^gauge T ^mediators T (3) 

where 

£gauge = -\g» v - \f )iv C x = ixH - m x xx (4) 

and the field strength tensors for the gluon and photon are G^ v and i 7 ^, respectively. Also, depending on the effective 
scenario the mediator lagrangian is 

^mediators = (Df i )*Df i -V i (5) 

where i = L,R and the covariant derivative, D, depends on the quantum numbers of /r or /r. The potential is 

Vi = \rn). f*fi + (Ai f* xPif + h.G.) (6) 

We remain agnostic about the mechanism that generates the x mass and assume it can be integrated out. Because 
we chose the Yukawa-type coupling between the dark matter, mediator and fermion to be non-zero, our analysis 
is generic and the derived bounds subsequently weaker in comparison to [421 which considers only the QCD gauge 
coupling between those particles. Once the quantum numbers of the sfermion are chosen, our scenario then has three 
basic independent parameters, 


{A,m/ i ,m x }. (7) 

These parameters are constrained by direct detection experiments, collider searches (see the relevant direct detection 
and collider citations in the Introduction) and unitarity constraints. Since we assume thermal dark matter and 
require x t° comprise some or all of the observed DM in the universe, our effective models are also constrained by 
the DM relic abundance. 


We consider scenarios with mediators that are color singlets as well as color triplets. The mediators may also 
be electroweak doublets as well. Although the perturbative unitarity bounds do not depend on the sfermion quantum 
numbers, the effects of these quantum numbers on the relic density can be reproduced through the following rescaling 

xrxPif^N^xrxPif ( 8 ) 


where the rescaling factor N depends on the representations of the groups fi is charged under. More details about how 
to compute N are given in Sec. |IV| Here N encodes the dimensions of the representations of the gauge group under 
which the mediator and SM fermions are charged. As N is increased, in general larger dark matter and mediator 
masses are possible because the dark matter annihilation is more efficient. Thus, there are weaker constraints from 
the relic abundance. As we show in Section [TV| our basic bounds scale accordingly. We briefly describe each scenario 
below. 


A. Colored Flavor Triplet Effective Model 


The colored mediators can have the quantum numbers 

(3,2)i/e (8; 1)— 2/3 (3, l)i/ 3 - (9) 

We consider all three possible types of mediators separately and will call them gT, ur and cIr respectively. In order 
to prevent flavor violation, we will consider the case in which these mediators are degenerate flavor triplets, 



ur = 




Ql = 


( 10 ) 
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with the DM-sfermion-fermion coupling matrix A proportional to the identity matrix. The collider, relic and direct 
detection constraints (but not the perturbative unitarity constraints) have been previously studied in [41]. 

To make our work more comprehensive, in the Appendix we modify the above scenario to allow for non¬ 
trivial couplings between the squark mediators and the SM Higgs. These couplings are analogous to the A-terms in 
SUSY models. Thus after electroweak symmetry breaking (EWSB), left-right squark mixing is generated. We focus 
on a model with one left-handed squark doublet and one right-handed up-type squark Ur, with the right-handed 
down-type squark cIr decoupled. Of note, this scenario has also been studied in the context of SUSY for the 
squark-mediated annihilation of bino-like dark matter [43]. 


B. Color Singlet Flavor Triplets Effective Model 


In this effective scenario, the color singlet mediators must have the quantum numbers 


(1; 2)-1/2 (M)l. (11) 

We denote these mediators as Ir and &r and call them sleptons. Once again, we consider them in degenerate flavor 
triplets, 



with A again proportional to the identity. The annihilation cross-section for color singlet mediators is decreased 
with respect to the colored mediator cross section. Thus, the relic density will play a larger role in constraining this 
scenario. Direct detection experiments will not be able to constrain this scenario at tree-level. 


C. Stop Effective Model 

Natural SUSY requires, at minimum, a light stop to cancel the SM Higgs mass quadratic divergences generated by 
SM top loops [38] . Thus far, the LHC has not discovered the stop, and is setting increasingly stringent constraints 
on the stop mass mm- Thus, because of its significance to natural theories of new physics, we also consider the 
case where the right-handed (RH) stop, tR , is much lighter than the other squark flavors. Although we only focus 
on light RH stop, models where only the left-handed (LH) stops is light should be constrained in a similar way. The 
LHC and, if necessary, a future hadronic collider is essential to discover/constrain this scenario. 


III. UNITARITY CONSTRAINTS 

We consider an effective theory where a charged/colored scalar (sfermion) mediates dark matter annihilation. The 
only couplings allowed to vary freely are the DM-fermion-sfermion Yukawa couplings, \ where i — L^R from 
equation ([6|. These couplings drive all the DM scattering and annihilation processes and are directly constrained by 
perturbative unitarity. As discussed in the Introduction, there has been a long understood conceptual dependence of 
the WIMP dark matter paradigm on perturbativity. 

It is well known [47U49] that any theory known to all orders in perturbation theory gives amplitudes which 
lie on the Argand circle. However most perturbative unitarity computations are done at tree-level and the loop 
corrections push the value of the scattering amplitudes toward the Argand circle. If the scattering amplitudes of 
some processes in the theory are too large at tree-level, unitarity cannot be restored without violating perturbativity. 
[ Ml 1481 [49] have derived a simple geometric argument that allows the translation of a perturbativity bound into an 
upper bound on the tree-level scattering amplitudes in any theory. This upper bound applies to the eigenvalues of 
the partial wave components of the transition operator T. This operator is defined such that 


S= 1 + iT 


(13) 
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where S is the S-matrix. In this work, we require that the loop contributions to the scattering amplitudes amount to 
less than 40% of the tree-level values. According to 05105], this amounts to enforcing 



(14) 


where 7/- is the j-th partial wave component of T. Since the lowest modes generally give the strongest constraints, 
we only focus on the s-wave component of T (j = 0). 


In the following, we derive unitarity bounds for the scattering of the yy and ff pairs where / is any of the 
fermions associated to the scalars decribed in Section (nj We compute the bounds in the large c.o.m. energy limit 

5 » M x , Mp v , (15) 

where v is the electroweak vev. XX XX scattering is forbidden at tree level in our model. At the center of mass 
energies considered here, both the electroweak and strong couplings can be neglected compared to A. The amplitude 
associated to //—)►// can therefore be neglected. We are thus left with only one process to consider, 

XX -> ff, (16) 

which is entirely determined by Xl and/or Xr. If the DM is a Dirac particle, XX ^ ff can only occur in the t-channel, 

through the diagram shown in Figure [l] Using the results quoted in f5Q| l5T ] . the amplitude for this diagram in the 

large s limit is 


M. =—X ^ 2 5(a x ,af)5(apaj) (17) 

where A^ is either A^ or A#, 5 is the Kronecker symbol and otf is the sign of the helicity of a given fermion /. Here, 
in the large s limit, we denote the helicity of the LH fermions as — and the one of the RH fermions as +. Since we 
work at large s, we can operate in the interaction eigenstate basis. For generic sfermions /#, /l, we then consider the 
following fermion pairs 


XlXl , XrXRi f lJl, f r$r. 


(18) 


and find 


T° = —— 
167r 


/ 0 0 0 X%\ 

0 o A| 0 

0 A| 0 0 

\An 000 


(19) 


The eigenvalues of this matrix are ±A^ and ±A|. For 0 G [0, j], the unitarity bound equation (14) then becomes 

|Al| , |Aft| < V8n ~ 5. (20) 


In the case of Majorana DM, both the t and u-channel diagrams contribute to the \X ff process (see Figure [l]). 
These diagrams interfere constructively and the resulting amplitude is doubled with respect to the Dirac case. The 
bounds quoted above then become 


|Al| , |A_r| < V47T ~ 3.5. 


( 21 ) 


IV. DARK MATTER CONSTRAINTS 

The DM constraints are effectively determined by the coupling A i associated to the DM-fermion-sfermion vertex. 
See equation §). This vertex parameterizes DM annihilation and, for light flavor squarks, DM-nucleus scattering at 
tree-level. 
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FIG. 1: Diagrams contributing to the XX ff amplitude for Dirac and Majorana DM. For Dirac DM, only the t-channel 
(top) diagram contributes to the scattering amplitude. 


A. Dark Matter Relic Abundance 


Dark Matter can annihilate to fermions via the diagrams shown in Figure [l] for Dirac and Majorana DM. For Dirac 
DM, annihilation occurs dominantly in the s-wave, with the s-wave term a being 


a oc N c Nf{ 

avors 


A 4 

87 t ( ra | - m 2 -) 2 


( 22 ) 


where Nc is the color factor. When the DM and sfermion masses increase, the coupling A also needs to increase to 
ensure efficient DM annihilation, unless the DM and sfermion masses are nearly degenerate. The coannihilation region, 
where the DM and mediator masses are within a few percent of each other, is however extremely narrow and can only 
be obtained at the price of a significant fine-tuning. Throughout this paper, we will define this coannihilation region 
as the region where the DM and sfermion masses are within at most 5% of each other. Outside the coannihilation 
region, the s-wave component of the DM annihilation cross section is bounded by 


A 4 

^max — 100 Nc N flavors —V • (23) 

m x 

In a perturbative theory, unitarity sets an upper bound on A, as studied in Section IH When A reaches its maximum 
value, a max will decrease as the Dark Matter mass increases. The unitarity constraints on A therefore translate into 
an upper bound on m x . This upper bound scales as y/NcNfi avors , independently on the other properties of the 
mediators involved. A similar reasoning leads to a neighboring upper bound on the mediator masses. 


A similar mechanism is at play for Majorana DM. In this case, however, DM annihilates through both the t 
and the ^-channels, as shown in Figure [l] The two corresponding diagrams interfere destructively. The s-wave 
component of the DM annihilation cross section to two fermions / is therefore suppressed by a factor of mj. Heavy 
DM annihilation then occurs mostly in the p-wave [52] . As computed in [53], the annihilation cross section for 
Majorana DM goes as 


^"ann 


E a4 


m< 


„2 m l 


(24) 


In the early universe, the velocity suppressed p-wave term dominates and, for degenerate sfermion masses, all the DM 
annihilation channels contribute equally to the relic density. This cross section is in general much smaller than the 
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FIG. 2: Example of coannihilation diagrams for the dark matter x and the sfermions /. Here, the wavy lines represent any 
gauge boson, including the gluon. 


annihilation cross section for Dirac DM. Outside the coannihilation region, the relic density constraints on Majorana 
DM are then expected to be significantly tighter than the constraints on Dirac DM. 

Inside the coannihilation region, Dark Matter can coannihilate with the sfermions through the diagrams shown in 
Figure [ 2 ] f54H57] . This process can significantly decrease the DM relic density and will considerably loosen the bounds 
obtained through the method described above. Note that the Griest and Kamionkowski bound [I] does not take 
coannihilation into account so, in this region, DM can be much heavier than 120 TeY without overdosing. 


B. Direct Detection 


DM-nucleus scattering occurs in the 5 and t (and u for Majorana DM) channels through intermediate squarks. In this 
work, we consider only tree-level DM-nucleus scattering and therefore focus on the light flavor squark components of 
ur cIr and gz,. The spin-independent (SI) direct detection (DD) cross sections for Dirac DM for ur (Ir and gx have 
been computed in m and are approximately 


a 


UR 

SI 


a 


dR 

SI 


r \j 


H 2 A 4 
64-7T (to 2 — to? ) 2 

v X UR / 

A 4 

64-7T (to 2 — TO? ) 2 

* OR 



Al „ 9m 2 A 4 
SI 647r (m| — m? t ) 2 


(25) 

(26) 
(27) 


Constraints from DD are in tension with relic density constraints since they prevent A from being too large. 

For Majorana DM, the operator responsible for spin-independent DM-nucleon scattering is generated at higher order 









in the effective operator expansion G2H5EH5D]. Outside the coannihilation region, the SI cross section is of the form 


v si ~ 


4 2 

m^m x 


(28) 


The SI direct detection cross section for Majorana DM is therefore strongly suppressed for most of the parameter 
space and will not lead to significant bounds. Although the spin-dependent (SD) cross section is not suppressed, 
the current bounds from SD direct detection experiments m are too loose to lead to significant constraints on the 
DM and mediator masses. For models involving both left and right-handed squarks, introducing a left-right mixing 
generates a dimension six effective operator contributing to the spin-independent DD cross section. Obtaining this 
left-right mixing in a natural way without inducing flavor mixing, however, leads to mixing angles that are strongly 
suppressed for light squark flavors. We show in the Appendix that introducing these small mixing angles does not 
significantly increase the spin-independent DD cross sections for models with Majorana DM. 


V. RESULTS 

In this section, we give the results when DM couples to the mediators: f ={u R , d R , qL , Ll, Lr}. For each case, we 
perform a uniform scan over 

Xi e [0, 5] m x , m~ u e [0,100] TeV (29) 

where i = L,R. We compute the DM relic density and direct detection cross section using MicrOmegas [62] and 
generate the model files using FeynRules 63 interfaced with CalcHEP [64]. We require the relic density to be less 
than the value measured by the Planck satellite [5] plus 3 cr, i.e., 

h 2 n c < 0.1226. (30) 

We require the points in the model to pass the LUX bounds |28]. We also check whether or not the surviving points 
are within the projected reach of XENON1T [26]. 


A. Results for Colored Flavor Triplet Effective Model 


Figure |3jshows the points that pass the unitarity, relic density and direct detection constraints in the (m x ,raj.) space 

for ur , d R and qL with Dirac DM. The narrow region along the diagonal is the coannihilation region, in which the 
relic density bounds are considerably loosened. Outside this region, the bounds on the DM and squark masses are 
around 60 TeV for u R and dz, and around 80 TeV for qL for which there exist s ix a nnihilation channels instead of 
three. Note that these upper bounds vary as V^flavors W, as explained in Section |rv[ Figure [I] shows the points that 
pass the unitarity, relic density and direct detection constraints in the (m x ,mj) space for u R (top plot), d R (middle 
plot) and qL (lower plot) with Majorana DM. As expected, the relic density constraint leads to tighter bounds than in 
the Dirac case. The bounds outside the coannihilation region are now around 10 TeV for all squark families. Figure [5] 
shows the direct detection cross section vs the DM mass for the points that pass the unitarity, relic density and the 
LUX direct detection constraints. For each parameter point, we rescale the direct detection cross section by the ratio 
of the relic density of x over the total DM relic density shown in equation (30). The points in red represent the points 


not generated by co-annihilation that are outside the projected reach of XENON1T. The yellow points are solely in 
the co-annihilation region. There were select by requiring 


m x — vrij 


< 5 %. 


(31) 


Direct detection experiments are extremely sensitive to the squark parameter space for dirac DM and the future 
experiments will be able to exclude almost all of our models, even in the case of co-annihilation. 


We also consider co-annihilating points in parameter space that satisfy equation (31). These points can sur¬ 
pass the well known bound from Griest and Kamionkowski [1] since it is not valid in case of coannihilation. Figure [6] 
shows the points that pass the unitarity, relic density and direct detection constraints in the (A, m x ) plane for u R , 
d R and qL- Both Dirac and Majorana DM are considered. We find that these bounds are around 300 TeV for the 
squarks. 
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FIG. 3: Points that pass the unitarity, relic density and direct detection constraints in the (rn x , my) space for uit (top plot), 
d,R (middle plot) and (]l (bottom plot) with Dirac DM. The narrow region along the diagonal is the coannihilation region. The 
colors denote the magnitude of the coupling of A i. See equations © and ( [29] ) . (Color online.) 
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FIG. 4: Points that pass the unitarity, relic density and direct detection constraints in the (m x ,m 
dR (middle plot) and qL (bottom plot) with Majorana DM. The narrow region along the diagonal is 
The colors denote the magnitude of the coupling of A i. See equations § and ( [29] ) . (Color online.) 


^) space for ur (top plot), 
5 the coannihilation region. 
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FIG. 5: Spin-independent direct detection cross section vs m x for the points that pass the unitarity, relic density and direct 
detection constraints for ur , cIr (top left and right) and qL (bottom) with Dirac DM. For each parameter point, we rescale 
the direct detection cross section by the ratio of the relic density of x over the total DM relic density shown in equation (30). 


The yellow points are the points in the coannihilation region while the red point are the points that are outside the reach of 
XENON1T. See equation ([3l|. The black dashed line represents the projected XENON1T bound extrapolated for large DM 







12 




2,5 0 100 200 300 400 500 

m x , m^ R , (TeV) 




m x , m~ dR (TeV) m x , m~ dR (TeV) 

FIG. 6: Points in the coannihilation region that pass the unitarity, relic and DD constraints for q^ (top), cIr (middle), and ur 
(bottom) for Dirac (left) and Majorana (right) DM. The horizontal line is the unitarity bound on the coupling A. 


B. Results for Stop Effective Model 

We plot the points that pass the unitarity, relic density and direct detection constraints in the (m xl m^ R ) space for 
RH stops, for Dirac and Majorana DM in Appendix [Bj Figure [9j For this model, we find bounds of 35 TeV for Dirac 
DM and 5 TeV for Majorana DM. Since we consider only one flavor here, the bounds on the DM and stop masses 
are tighter than the squark bounds by a factor of about \/3. Since the top component in the nucleons is largely 
subdominant, no significant DD bound can be obtained. 
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FIG. 7: Important squark pair production modes. For the large coupling A considered here, the t-channel diagram at 
the bottom of the figure can contribute significantly to the squark pair production cross section. Thus, if these models are 
responsible for the measured relic abundance, the sensitivity of a 100 TeV hadron collider can be greatly enhanced. 


C. Results for Colored Singlet Flavor Triplet Effective Model 


We plot the points that pass the unitarity, relic density and direct detection constraints in the (ra x ,raj) space for 

e,R and U with Dirac and Majorana DM in Appendix [bJ Here the bounds are tighter than for the squarks since the 
slepton annihilation is not enhanced by the color factor. The mass bounds are therefore decreased by a factor of y/3. 


For Dirac DM, we find a bound of 35 TeV for ep and 50 
TeV for e# and U, respectively. Also Appendix [Bj Figure 
which satisfy equation (31). 


eV for II. For Majorana DM, we find a bound 5 and 9 
has the slepton results for the co-annihiliating points 
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D. Comments on Collider Constraints 

Majorana DM models as well as all models with only stops or sleptons are not significantly constrained by direct 
detection. However, both DM and sfermions can be searched for at colliders. In SUSY, the dominant production 
mechanism for squarks is pair-production with an intermediate gluon. In our study, A is required by relic density 
constraints to be large for heavy DM, and t-channel squark production through the diagram shown in Figure [7] starts 
becoming significant. Bounds on the ur , (Ir and gz, models using the 8 TeV ATLAS and CMS data have been derived 
in [41] • These bounds are about a few hundreds of GeV and exclude only a tiny fraction of the parameter space 
considered here. A 100 TeV hadron collider, however, could exclude squarks up to a few TeV, as shown in 1551 166 ] . 
Such a collider could exclude a sizeable part of the parameter space, especially for Majorana DM. For DM heavier 
than a few TeV, the relic density requirement constrains the masses of the DM and the squarks to be close to each 
other, and traditional squark searches start losing their power. Due to the size of the coupling A, mono jet searches 
are a serious competitor to squark searches in these kinematically challenging regions. 


VI. CONCLUSIONS 


Understanding the nature of dark matter is one of the most compelling question in particle physics. Dark matter 
that was once in thermal equilibrium with the SM is very well motivated. In this work we focus on scenarios 
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where the thermal dark matter annihilates via new charged and/or colored scalars (sfermions). We have shown 
here that combining unitarity arguments with relic density constraints allows to derive upper bounds on both the 
masses of the DM and of the sfermion mediators. Using flavor arguments, we consider various possible sfermion 
configurations, notably, degenerate sfermion flavor triplets and models where the stop is much lighter than the other 
new particles. For Dirac DM, the upper bounds on the DM and sfermion masses vary between about 35 TeV for stops 
and right-handed sleptons to 80 TeV for left-handed squark doublets outside the coannihilation region. The bounds 
for Majorana DM are much tighter, varying from 5 TeV for stops and right-handed sleptons to about 10 TeV for 
left-handed squark doublets. We also consider co-annihilation regions, in which the Griest and Kamionkowski bound 
does not apply. In these regions, we derive upper bounds on the DM mass as high as 400 TeV. The upper bounds 
derived in this work strongly depend on the number of scalar mediators driving the DM annihilation. This work, 
however, spans a wide range of possible configurations and, within the MFV framework, can be readily adapted to 
models involving different combinations of sfermions. 

Exploiting the tension between unitarity and relic density constraints on DM and sfermions not only allows 
us to derive upper bounds on the scales of new physics, but also to define goals for future experiments. For Dirac 
DM, in particular, future direct detection experiments such as XENON1T could probe most of the parameter space 
in models with light-flavor squarks, leaving only extremely narrow regions viable. For Majorana DM or models with 
stops or sleptons, however, collider searches and indirect detection will provide the best constraints on new physics. 
For Majorana DM, which is associated with particularly tight bounds, a future lOOTeV collider would be able to 
probe most of the parameter space for squark models. Since relic density constraints require the DM and sfermion 
masses to be close to each other, however, developing and improving compressed spectrum searches will be crucial to 
get a better understanding of scalar mediator models. 
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Appendix A: Squark Flavor Triplets with Left-Right Mixing 

If we allow for a non-trivial coupling between the sfermion mediators and the SM Higgs that is analogous to the 
A-terms in SUSY, we can generate left-right (L-R) mixing after electroweak symmetry breaking (EWSB). This can 
be important for models with flavor triplet squarks and Majorana DM. In this case, the spin-independent direct 
detection cross section is suppressed due to its twist-2 origin. The direct detection constraints could in principle 
become stronger if the left and right-handed squarks are allowed to mix. This scenario has been most recently studied 
in the context of SUSY for sneutrino-mediated DM annihilation [33]. To explore this possibility, we consider a model 
with one left-handed squark doublet gz, and one right-handed up-type squark Ur, with the right-handed down-type 
squark <1r decoupled. Since mixing between different squark flavors is strongly constrained by kaon oscillations m, 
we only allow mixing between left and right-handed squarks of the same flavor. Such a scenario can be obtained 
naturally by requiring the minimum flavor violation (MFV) [44jf45], an d requiring the left-right squark mixing terms 
in the Lagrangian to be proportional to the masses of the corresponding quarks. The mass and interaction terms for 
the squarks are then taken to be 


Vmass = \rnl L q* L q L + ^m\ R u* R u R 

Vint = ^LQlXPlQ + ^RUrXPrU + h.c. 
VVukawa = A LR Y u q* L HuR + h.c. 


(Al) 

(A2) 

(A3) 


where Y u is the matrix of the SM up-type Yukawa couplings, H is the SM Higgs, Alr is an A-term-like parameter 
with dimensions of ene rgy, and the squark mass matrices rriq L and rriu R are proportional to the identity. The Yukawa 
potential in equation (A3) induces a mass mixing between the up-type component of gz, and ur after EWSB. The 
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m x (TeV) 


FIG. 8: crsi vs m x for points with 63 < 0.1 (blue) and O3 > f —0.1 (green). The points displayed here have passed the unitarity, 
relic density and LUX constraints. 


mass eigenstates u\ and fq are then such that 

ul = cos# ft 1 + sin0u2 (A4) 

ur = — sin# ft 1 + cos0u2- (A5) 

Since the L-R mixing term from the squark-Higgs interaction is proportional to the SM Yukawa couplings, the mixing 

angle 6 is flavor-dependent, with the stops mixing much more strongly than the charm and up squarks. Since we have 

decoupled d#, the LH down squarks do not mix with anything. The DM interaction terms with fq, U 2 and dj, are 
then 


Vint = u*i X (Al COS 9 P l - \r sin 9P R ) u + u * 2 \ (a l sin 9 P L + X R cos 0 Pr) u + d L \PLd + h.c. (A6) 

Unlike the previous effective scenarios, the resulting effective model has six free parameters, 


{m x ,rriq L ,mu R , A LR , \ L , X R } . (A7) 

For Majorana DM, introducing left-right mixing generates dimension six effective operators that contribute to the 
spin-independent Dark Matter direct detection cross section. For a given squark flavor i, the lowest order mixing 
contribution to the SI cross section is of the form [43] 


c 7 si oc X 2 l X\ sin 2 2 9 


mi 


mi 


(AS) 


Due to the hierarchy of the SM Yukawa couplings, sizeable mixing will occur only in the stop sector. For light flavor 
squarks, the enhancement of the DD cross section due to introducing left-right mixing is therefore expected to be 
negligible. Fig.[8]shows the points that pass the unitarity, relic and LUX constraints in the (ra x , <tsi) plane. The blue 
points are points for which the stop left-right mixing is low 


Os < 0.1 

while the yellow points are points for which the stop left-right mixing is large 

0s> J-0.1. 


(A9) 


(A10) 
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FIG. 9: Top plots: Points that pass the unitarity, relic density and direct detection constraints in the pace for 

t R with Dirac (left) and Majorana (right) DM. Middle plots: Points that pass the unitarity, relic density and direct detection 
constraints in the (m x ,mj) space for e R (left) and II (right) with Dirac DM. Bottom plots: Points that pass the unitarity, 
relic density and direct detection constraints in the (m x ,mj) space for e R (left) and II (right) with Majorana DM. For each 
plot, the narrow regions along the diagonal is the co-annihilation region. 
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m x , m/ (TeV) m x ,mi L (TeV) 


FIG. 10: Points in the coannihilation region that pass the unitarity, relic and DD constraints for e,R (top) and II (bottom) for 
Dirac (left) and Majorana (right) DM. The horizontal line is the unitarity bound on the coupling A. 


We removed the points belonging to the coannihilation region by requiring 


m x - m^s 

m x 


> 10 %, 


(All) 


ra| being the mass of the lightest stop, which is also the lightest squark in our model. As expected, introducing left- 
right mixing does not lead to a significant global improvement of the spin-independent DD cross section for Majorana 


DM. 


Appendix B: Additional Results for Various Effective Models 

Considering the effective field theory with only stops and dark matter, we plots all of the points that pass the 
unitarity, relic density and direct detection constraints in Figure M The plots are in the (m x ,m^) space for RH 
stops, for Dirac and Majorana DM. Because there is only one squark flavor, the bounds on the DM and stop masses 
are tighter than the squark bounds by a factor of about y/3. 

Figure [9] also shows the unitarity, relic density and direct detection constraints for the effective field theory 
with flavor triplet sleptons. Again, the plot is in the (m x ,mj) space for e r and II with Dirac and Marjorana DM. 
Again, the bounds are tighter than for the squarks bounds shown in the text since the slepton annihilation is not 
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enhanced by the color factor. The mass bounds are therefore decreased by a factor of y/3. Since 7Vfl avor N c is the 
same for the e# model (reminder, three flavors) and the stop-only model, the bounds on the DM and squark masses 
are the same in both cases. 

Finally, we consider the parameter space points in the co-annhilation regime for the slepton effective theories. 
Figures [6] and [lO] shows the points that pass the unitarity, relic density and direct detection constraints in the 
(A, ra x ) plane for II and e#. We consider both Dirac and Majorana DM. As emphasized before, since Griest and 
Kamionkowski [l[ is not valid in case of coannihilation, these bounds can be well beyond 120 TeV. We find 200 TeV 
for RH sleptons for both Dirac and Majorana DM. 
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